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ABSTRACT 

We present redshifts and optical richness properties of 21 galaxy clusters uniformly selected by 
their Sunyaev-Zel'dovich signature. These clusters, plus an additional, unconfirmed candidate, were 
detected in a 178 deg 2 area surveyed by the South Pole Telescope in 2008. Using griz imaging 
from the Blanco Cosmology Survey and from pointed Magellan telescope observations, as well as 
spectroscopy using Magellan facilities, we confirm the existence of clustered red-sequence galaxies, 
report red-sequence photometric redshifts, present spectroscopic redshifts for a subsample, and derive 
R200 radii and M200 masses from optical richness. The clusters span redshifts from 0.15 to greater 
than 1, with a median redshift of 0.74; three clusters are estimated to be at z > 1. Redshifts inferred 
from mean red-sequence colors exhibit 2% RMS scatter in <j z /{1 + z) with respect to the spectroscopic 
subsample for z < 1. We show that M200 cluster masses derived from optical richness correlate with 
masses derived from South Pole Telescope data and agree with previously derived scaling relations to 
within the uncertainties. Optical and infrared imaging is an efficient means of cluster identification 
and redshift estimation in large Sunyaev-Zel'dovich surveys, and exploiting the same data for richness 
measurements, as we have done, will be useful for constraining cluster masses and radii for large 
samples in cosmological analysis. 

Subject headings: galaxies: clusters: individual, cosmology: observations 
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1. INTRODUCTION 

Galaxy cluster s are laborato ries for both astrophysics 
and cosmology ( |Evrard||2004j ). Clusters represent the 
most massive dark matter halos, and their number den- 
sity as a fun ction of cosmic time is highly sensitive to 
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dark energy (Wa ng fc Steinhardt||1998| |Haiman et al. 
2001] [Holder et al.| 1^0011 l^attye & Wellerl 1^00^1 |Mol- 
nar et al.||2004| |Wang et al ||2004| |Lima Hu||2007p . 
The mass ol these systems is dominated by dark matter, 
but the primary means of observing clusters — especially 
large samples of them — are the luminous baryons of the 
hot intracluster gas and the galaxies themselves. The 
formation of the halos is well understood, while the pre- 
cise behav ior of the baryons is not as well modeled (see 
|Voit|2QQ5| for a review). This gap must be closed so that 
data from large cluster surveys can place precise con- 
straints on cosmological parameters over a wide range 
of redshifts. Mult i- wavelength observations of a cleanly 
selected, redshift-independent sample of galaxy clusters 
are a potentially powerful method of achieving this. 
Searches for galaxy clusters using the Sunyaev- 



ZePdovich (SZ) effect flSunyaev fc Zerdovich| |1972 ) 
promise to provide such a clean, redshift-independent 
sample. The SZ effect is scattering of cosmic microwave 
background photons to higher energy by the hot elec- 
trons in galaxy clusters ( |Birkinshaw||1999| ) . The SZ sur- 
face brightness is independent ol redshift but is closely 
related to cluster mass and so it is expected to be an ex- 
cellent method for creating approximately m ass-limited 
samples ext ending over a wide redshift range (|Carlstro"m 



et al. 2002). The constraints on cosmological parame- 
ters from such samples are complementary to geometri- 
cal tests using ty pe la supernovae and baryon acoustic 
oscillations (e.g., Vikhlinin et al. 2009| ). Two SZ sur 
veys, the A tacama Cosmology Telescope (ACT; | Fowler 



et al.||2007[) and the South Pole Telescope (SPT; |Carl- 



strom et al.j2009 ) projects, are well positioned to provide 



large surveys which can be used for growth of structure 
studies. 



Staniszewski et al. (2009, hereafter S09) presented the 
first discovery ol previous unknown galaxy clusters using 
their SZ signature. Cluster redshifts are needed in ad- 
dition to SZ data to provide the strongest constraints 
on dark energy. Coordinated optical follow-up obser- 
vations can provide the needed redshift measurements 
The Blanco Cosmology Survey (BCS; Ngeow et al 
2006 and http : / / cosmology . Illinois . edu/BCS/ ), ar 



IN (J AO survey program (2005-2008), provided multiband 
optical observations for the initial follow-up of portions of 
the first SPT survey fields. These data were used to iden- 
tify optical counterparts to the |S09| sample, search for 
giant arcs, explore possible cluster superpositions, and 
derive photometric redshifts. 

Cluster mass can be estimated using several methods: 
the SZ and X-ray luminosity, which are sensitive to in- 
tracluster electrons; the number, luminosity, and veloc- 
ity dispersion of cluster galaxies; and from gravitational 
lensing , which is the most dir e ct pro be of total cluster 



Menanteau & Hughes ( [2009) ch aracterized the 
gala xy counts and luminos ity of the S09 cluster sample, 
and |McInnes et al.| ( |2009[ ) subsequently explored their 
weak gravitational lensing signals. 



et al. (2010 hereafter V10) present an additional 17 SZ- 
detected clusters. Here we describe coordinated optical 
imaging of the catalog of 21 uniformly selected SZ detec- 
tions, and new spectroscopic results on 8 of the clusters. 
Counterparts to a subset have been found in t he cata- 
logs of [Abell et al | ( |1989[ hereafter A89) and |Menan-| 
teau et al.| (|2010[ hereafter SCS-II). Seven clusters tell 
within the BCS footprint. For the remaining 14 clusters, 
and also for a subset of the BCS sample, we conducted 
pointed imaging observations and, for 8 clusters, spectro- 
scopic observations, with the Magellan telescopes. The 
photometry was used to search for overdensities of red- 
sequence galaxies near the SZ locations, and if present, 
estimate their redshifts and also characterize their mass 
via optical red- sequence galaxy counts, or richness. 

We describe in £[2] the observations and data reduction. 
Section [3] outlines the redshift and richness an alys is we 
used, and jp] d escribes the results on redshift (£4.1) and 
richness (j |j~2| ). In Section [5] we discuss the results, and 
conclude with ^6] 

Throughout this paper we assume a flat concordance 
ACDM universe with (O a ,O m ,^) = (0.736,0.264,0.71) 
(Dunkley et al. 2009). All magnitudes are in the Sloan 
Digital Sky Survey (SDSS) griz AB system. 

2. DATA ACQUISITION & REDUCTION 

Cluster detection was achieved using millimeter- 
wavelength data from the South Pole Telescope, and op- 
tical imaging and spectroscopy provided cluster confir- 
mation and redshift and richness estimates. 



2.1. South Pole Telescope 



The sample of 21 clusters presented here and in V10 
is the first cosmologically significant catalog of clus- 
ters selected via the SZ effect. The sample was se- 
lected from two SPT survey fields totaling 178 deg 2 at 
R.A. = 23 h 30 m , Decl. = -55° and R.A. = 5 h 30 m , 
Decl. = -55° (J2000). Both fields were observed with 
arcminute resolution to an equivalent white noise level 
of 18 /iK-arcmirpI 

The SPT time-ordered data were filtered and binned 
into maps, with the filtering acting roughly as a 1° high- 
pass filter in the R.A. direction. Clusters were extracted 
from these maps using a matched filter approach based 
on the work o f (Haehnelt fc T egmarkj |1996| Herranz 
et al.||2002a|b| |Mehn et al.||2006p . SpakiaTTilters were 



Using data acquired by the SPT in 2008, Vanderlinde 



constructed to maximize detection significance within a 
set of cluster profiles. The SPT astrometry is based on 
comparisons of radio source positions derived from SPT 
maps and positions of tho se sources in the AT20G cat- 
alog (Murphy et al. 20101), and should be accurate to 

5". 

Cluster candidates were then identified by selecting all 
peaks above a fixed significance threshold and choosing 
the filter scale which produced the maximum detection 
significance £. A3 parameter model fo r M20 involving 
£ and the cluster redshift is presented in |V10[ along with 
the details of the SPT data reduction. 

31 The unit K refers to equivalent fluctuations in the CMB tem- 
perature, i.e., the level of temperature fluctuation of a 2.73 K black- 
body tha t wou ld be required to produce the same power fluctua- 
tion. SeelVTo] 
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2.2. Blanco Cosmology Survey 

BCS is an NO AO survey program to obtain deep 
griz imaging of two southern fields centered at R.A. = 
23 h 00 m , Decl. = -55°12 / and R.A. = 5 h 30 m , Decl. = 
-52°47 / (J2000), each roughly 50 deg 2 . The 2008 SPT 
survey fields are larger than the BCS fields and include 
the entire BCS regions. BCS was conducted from 2005- 
2008 using the Mosaic-II wide-field imager on the Blanco 
4 m telescope at Cerro Tololo Inter- American Observa- 
tory. Mosaic-II is an array of eight 2 k x 4k CCDs with a 
pixel scale of 0.270" pixel -1 and a 0.36 deg 2 field of view. 
The strategy was to obtain deep, contiguous griz imag- 
ing of the survey fields. In addition, BCS imaging was 
carried out for 7 photometric redshift calibration fields 
which include a sample of several thousand published 
spectroscopic redshifts. 

The BCS data were processed and calibrated using a 
data manage ment system developed for the Dark En- 
ergy Survey (Ngeow et al. 2006| |Mohr et al~ 2008) and 
run on the NCSA TeraGrid !A-o4 Linux cluster. Data 
reduction includes crosstalk correction, overscan correc- 
tion, bias subtraction, flat fielding, fringe and illumina- 
tion corrections, field distortion correction, standard star 
photometric calibration, coadd-stacking, and photomet- 
ric extraction of sources. BCS stacks typically reach 5cr 
galaxy photometry limits of 24.75, 24.65, 24.35 and 23.5 
mag in griz, corresponding, for example, to a 0.5L* clus- 
ter elliptical at z = 1. Completeness measurements in 
a typical field that are derived from comparison with 
deeper, better seeing CFHT data suggest 50% complete- 
ness limits of 24.25, 24.0, 23.75 and 23.0 in griz (Zenteno 
et al. in prep.). 

2.3. Magellan 
2.3.1. Imaging 

For clusters that fell outside the BCS coverage re- 
gion, for 5 that were within the BC S region, and for 



the unconfirmed candidate (see [4.3) we obtained griz 
imaging with the Inam ori Magellan Areal Camera and 
Spect rograph (IMACS; iDressler et al. 1120031 lOsip et al. 



2008 ) and Low Disper sion Survey Spectrograph (LDSS3; 
see TDsip et all 120081)— both in imaging mode — on the 
twin Magellan 6.5m telescopes. IMACS is on the Mag- 
ellan Baade telescope at the //ll Nasmyth focus. Its 
circular field of view in f/2 mode subtends 0.20 deg 2 , 
mapped onto a 8192 x 8192 pixel, 8-chip CCD array, for 
a pixel scale of 0.200" pixel -1 . LDSS3 is on the Mag- 
ellan Clay telescope at the //ll Nasmyth focus, and its 
roughly 60 arcmin 2 circular field of view maps onto a sub- 
region of a 4064 x 4064 pixel CCD, at a pixel scale of 
0.189" pixel -1 . 
Using the sci ence-tested pipeline de scribed by |Rest 



et al, 



(2005) and 



Miknaitis et al. 



(2007), the same image 



reduction operations described above were performed on 
the Magellan images. Instead of using standard stars for 
photometric calibrati on, however, we e mploy Stellar Lo- 
cus Regression (SLR; High et al. 2009| ) using stars that 
appear in the cluster images themselves. SLR delivers 
photometric calibrations by regressing the instrumental 
color-color locus of stars in any field to the known, as- 
trophysically fundamental locus in the AB system. SLR 
enabled us to forego observations of standard stars al- 



together, maximizing the total integration time on the 
cluster fields. 

The principal difference in strategy between Magel- 
lan and BCS observations is that the Magellan exposure 
times were adaptive rather than uniform. We first ex- 
posed for roughly 100 s in griz, searched for a cluster 
in the images, and continued with additional exposures 
if none was found. In median seeing of 0.8", these ex- 
posures reach nominal, 5<j point-source limiting magni- 
tudes of 24.8, 24.8, 24.4, and 23.4 mag in griz. If no 
cluster was detected with the initial set of images, we ac- 
quired further exposures until a detection was achieved 
at a depth of approximately 1 to 0.4L* with respect to 
the early- type cluster galaxies. This strategy results in 
highly variable depth for the Magellan imaging, but is 
the most efficient use of telescope time for follow-up ob- 
servations. 

Combined optical and SZ images are shown in Ap- 
pendix |Aj Figures [7}|28l 



2.3.2. Spectroscopy 

Spectroscopic data were acquired with LDSS3 in long- 
slit mode for the purpose of measuring redshifts. Given 
limited telescope time for spectroscopy, we observed only 
a subset of the confirmed clusters. The subset was chosen 
to span the widest possible range in redshifts, so that we 
could assess the performance of o ur re d- se quen ce redshift 
measurement methodology (see £3.2 and j |4.1| ). 

We obtained low-resolution spectra of galaxies in the 
field of 8 SPT clusters, using the VPH-Red and VPH-A11 
grisms. The median seeing was about 0.7" and condi- 
tions were photometric. Standard CCD processing and 
2D-spectrum extraction, with preliminary wavelength so- 
lutions, were accomplished with the COSMOS reduction 
packag^j the ID spectra were then extracted using the 
apall task in IRAF. We employed our own IDL routines 
to flux calibrate the data and remove telluric absorption 
using the well- exposed continua of the spectroph otomet- 
ric standards ( |Wade fc Ho"rne|[T988l [Foley et al.(2QQ3| |. 



3. ANALYSIS 

We adopted a standard red-sequence model for optical 
cluster detection, photometric redshifts, and richness es- 
timation. As previously introduced, spectroscopic data 
on about half the clusters were used to empirically cor- 
rect the red-sequence model colors and then verify pho- 
tometric red-sequence redshifts over a long redshift base- 
line. 



3.1. Red- Sequence Model 
Red-sequence models w ere derived from the work of 



Bruzual fc Chariot] (2003), reflecting passively evolving 
instantaneous- burst stel lar populations with a formation 
redshift of z = 3, using Berte lli et al. l ( 1994[ ) evolution- 
ary tracks and the |ChabrTe r (2003) initial mass function. 
At each redshift, a range of metallicities was chosen by 
including a randomization in the metallicity-luminosity 
relation. The models were smoothed by linear fits in 
color-magnitude space at each redshift, and finally inter- 
polated to arbitrary redshift using cubic splines. 

32 http:/ /obs. carnegiescience.edu/Code/cosmos/ 
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Fig. 1. — Red-sequence color-magnitude models as a function of 
redshift. Galaxy colors r — i and i — z were fit simultaneously in 
the red-sequence analysis, because that particular combination of 
colors (namely, r—z = (r—i)-\-(i—z)) allows for monotonic mapping 
from color to redshift over the widest redshift range (0 < z < 1.4) 
using optical wavelengths. The apparent i band magnitude of M* 
at each redshift is denoted with the black points. 

We made an empirical correction to the model using 
the subsample of 10 clusters with spectroscopic data . Af- 
ter performing the red sequence peak finding (£3.2), we 
plotted red sequence redshifts against spectroscopic red- 
shifts. The best fit line through the data was measured to 
have slope 0.89 ±0.03 and ^-intercept 0.04 ±0.02. Model 
colors were then corrected by reassigning the model red- 
shifts to be equal to the inverse of this linear relation. 
The effect of the model color correction is to leave red- 
shifts near the pivot, z ~ 0.4, roughly unchanged, and 
to boost redshifts near ~ 1 by about Az = ±0.07, or 
Az = ±3.5% (1 ± zY The resulting red-sequence model 
is shown in Figure [T] 

3.2. Red- Sequence Finding and Redshifts 

We confirm the existence of clusters by searching near 
the SZ cluster coordinate for a background-subtracted 
excess of red- sequence objects, effectively segregated by 
redshift. At each redshift step of size Az = 0.01 in the 
range 0.1 < z < 1.4, and with an aperture of 2' radius 
around the SZ coordinate, we select all objects with a 
photometric signal-to-noise ratio > 5 in riz, whose r — i 
and i — z colors are also within 2a of the red-sequence 
model line. The total uncertainty, e.g., for r—z, is defined 
as 

2 _ 2 , rs2 / 1 \ 

a = a ri + a ri , (1) 

and likewise for i — z. Here a r i (cr^) is the photo- 
metric uncertainty in r — i (i — z) color of an object, 
and a™ {a\ s z ) is the intrinsic color scatter of the red se- 
quence. These colors were chosen because their combi- 
nation, r — z = (r — i) ± (i — z), increases monotonically 
over a long baseline of redshift, < z < 1.4. We as- 



smaller than this, and has been shown to be constant 
with redshift out to z « 1.2 flMenci et al.|2008||Mei et al.| 
|2009| . Morphological and spectral galaxy classification 
is beyond the scope of this work, and we assume our 
red sequence selections contain SO galaxies in addition 
to ellipticals. This increases the color scatter and may 
also effectively induce redshift evolution due to evolving 
SO populations; we ignore the latter effect, as the small 
number of clusters presented here is not sufficient to con- 
strain redshift evolution. 

We sum the selected galaxies, and normalize the counts 
by the projected area. This yields total surface density of 
all objects with colors consistent with the red-sequence 
model, within the 2! aperture, as a function of redshift, 
Stotai(^)- Uncertainties are estimated as the square root 
of counts, divided by the area. 

The background red-sequence surface density is mea- 
sured in a similar way. In the same set of viz expo- 
sures, we count red- sequence objects within many ad- 
jacent apertures of size 5' x 5' over the entire field of 
view, excluding the 2! region around the SPT candi- 
date position. The background surface density at each 
redshift, ^backgrounds), is calculated as the median of 
area-normalized counts from all apertures, and the un- 
certainty is the standard deviation divided by the square 
root of the total number of apertures. This approach is 
meant to minimize the contamination of the background 
signal by the red galaxies in the clusters themselves. 

The red-sequence excess as a function of redshift is 

^net = ^total ^background (2) 

in units of galaxies per arcmin 2 . We then renormalize 
the red sequence excess S ne t in each redshift bin such 
that each galaxy contributes a constant area under the 
curve. The reason for this is that one galaxy may fall into 
multiple adjacent bins due to the size of color uncertain- 
ties, sometimes significantly widening the red-sequence 
peaks, especially for faint objects. The renormalization 
is a measure to mitigate this effect. 

A cluster is detected if (1) there is an excess of red 
galaxies of the same apparent color in false color im- 
ages, and (2) the galaxies corresponding to the maxi- 
mum overdensity, S net (z max ), are those identifiable in 
the false-color images. The cluster photometric redshift 
is taken to be £ max . Figure [2] illustrates the process of 
red-sequence finding for one of the SZ clusters. 

3.2.1. Completeness 

We tested for completeness in red-sequence finding at 
representative BCS and Magellan depths. After con- 
structing mock optical catalogs (Song et al. in prep.) 
containing simulated clusters of mass above ~ 3 x 
10 14 h~ x Mq, we searched for red-sequence overdensities 
as we have described. Our estimated completeness is 
the fraction of clusters recovered from the mock catalog. 
These mock catalogs include clusters with red and blue 
galaxy populations that are tuned to match the pop u - 
lations observed in real clusters (e.g., [Lin et al.| 2004). 
The galaxy distribution in space is determined using 
subhalo positions within high-resolution A^-body simu- 
lations, and so the mock cluster spatial, kinematic and 
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Fig. 2. — An illustration of red-sequence finding for the cluster SPT-CL J0551-5709 at z spe c = 0.4230. The top left plot shows the surface 
density of objects consistent within 2cr of the red-sequence model as a function of redshift. The peak occurs at redshift z m ax = 0.41, which 
is consistent with the spectroscopic redshift (vertical white line, with a z = 2%(1 + z) uncertainty region shaded in gray for illustration). 
The color-magnitude diagram for all objects within 2' of the SPT coordinate is shown in the upper right, and the subsample of red-sequence 
objects at z m ax shown as red points. The vertical dotted line is the model m* at this redshift. Positions of all objects consistent with 
the Zmax red sequence are shown in the bottom panel, where we have also circled the cluster aperture. A spatial overdensity of objects is 
clearly seen at the aperture. 

color signatures are a good match to those seen in real 
clusters. 

The resulting selection function is shown in Figure [3j 
We recovered 100% of simulated clusters above the given 
mass threshold up to redshift 0.9. At this point the 
4000 A break begins to redshift out of the i band, mak- 
ing galaxies much harder to detect in i — z color space. 
The completeness begins to fall here, and for these BCS 
and Magellan depths the probability of finding an opti- 
cal counterpart above this mass threshold falls to zero by 
z ~ 1.2. Near-infrared photometry or much deeper space 
based optical photometry is required to push reliably to 
higher redshifts. But for the current sample only a sin- 
gle SPT cluster candidate was not confirmed using this 
multiband optical method. 



3.3. TVgai and A200 Richness 

After confirming a cluster and estimating its redshift, 
we measured the optical richness using a p rocedure that 
emula tes the MaxBCG richness estimator (Koester et al. 
2007), but adapted to high-redshift clusters. This began 
by again selecting objects with color within 2a of the 
red-sequence line at redshift £ max ; luminosity brighter 
than 0.4L* and fainter than the brightest cluster galaxy 
(BCG); and position within a projected radius of R = 
1 h~ x Mpc of the cluster center. We binned the selected 
objects in i band magnitude bins of size Am = 0.4. We 
subtracted the background, as before, by performing the 
same procedure in many apertures on the sky away from 
the cluster but in the same exposures, and normalizing 
by projected area. The background was subtracted from 
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Fig. 3. — This shows the completeness of optical cluster finding 
from tests on a mock galaxy catalog with depth representative of 
the BCS survey (solid line) and Magellan imaging (dashed line). 
We have approximated Magellan data here as having limiting mag- 
nitudes ~ 1 magnitude brighter than BCS. 

the total red-sequence counts in each magnitude bin. 

We then fitt ed Schechter luminosity functions (LF, 
|Schechter|1976 ) to the i band magnitudes of the selected 
objects: 



i)dm = O.41n(lO)0*lO-°- 4(m - m * )(a+1) 



x exp 



_ ^ Q — . 4 ( rn — rn * ) 



dm, 



(3) 



where m* is the characteristic magnitude of the LF, a 
is the faint-end slope, and (jf is the normalization. Our 
photometry was not uniformly complete on all clusters 
with respect to m*, so our ability to constrain the slope 
was weak. We therefore fixed a = — 1, which has been 
shown to be reasonable for the most massive MaxBCG 
clusters (i Linfc Mohr 2003; Hanse n et al.|2005| |Rudnick 
et al.||2009[ |Crawlord et al.||2UU9| ). We tested fixing wT 
at each redshift using our passive model, and leaving it 
free. For those clusters for which the luminosity function 
was well constrained, we measured values consistent with 
the model, and we ultimately chose to fix this parameter 
at each redshift. This is in agreement with the detailed 
LF studies of Zenteno et al. (in prep) on a subset of 
these clusters. The f aint-end slope a has been shown to 
evolve with redshift QRudnick et al.|2009[) , however, such 
studies must be performed at magnitudes fainter than 
m* + 1. We test varying a by ±0.3, and find our richness 
results are largely unaffected to within our uncertainties. 
Because we integrate only to m* + 1 there is only weak 
sensitivity to the adopted faint-end slope. 

We fit the LF only at magnitudes brighter than our 
limiting magnitude and fainter than the BCG. Limiting 
magnitudes are estimated to be the faintest magnitude 
bin before which the red-sequence background becomes 
incomplete, as indicated by a deviation from linearity 
in the logarithmic i band magnitude distribution of red- 
sequence galaxies. We also assess photometric complete- 
ness using simulated point sources in our cluster images. 
We determine that we recover 90% of simulated objects 
in images at all magnitudes brighter than the above- 
defined limiting magnitude. The majority of unrecovered 
objects are lost to pixel masking due to bright stars and 



o.i 



R < 1 Mpc| 





0.1 



Fig. 4. — An illustration of luminosity function fitting for SPT- 
CL J0551-5709. We use our passively evolving model to fix m* 
to 19.77 mag (vertical dotted line), and integrate down to m* + 1 
to estimate N ga \ (top panel) and N200 (bottom panel). R200 is 
estimated from iV ga i, as described in the text. 

to object crowding. We make a generic correction by this 
amount in all magnitude bins, the effect of which is to in- 
crease 7V ga i estimates by 0.9 -1 . The analytic LF is then 
integrated down to m* + 1. The resulting integration is 
an estimate of iV ga i, which is the number of red sequence 
galaxies within 1 h~ x Mpc of the cluster center, above a 
fixed luminosity threshold. 

We stress that the LF fitting step deviates from the 
MaxBCG procedure, but it is necessary because our pho- 
tometry is not complete to 0.4L* on all clusters. Our 
fitting and extrapolating the LF of fairly bright, red- 
sequence satellite galaxies is physically reasonable be- 
cause it is known that the luminosity function of such 
members of massive clusters is only very weakly depen- 
dent on mass and is well described by a Sc hechter func- 
tion from lower redshifts ([H ansen et al. 2009) to redshifts 
close to unity (IGilbank et al.l|2008|) . Figure [4] illustrates 
the resulting LF fit lor one of the SPT clusters. 

R200 is estimated fr om 7V ga ] using the empirical N ga \- 
R200 relation of |Hansen et al.| ( |2005[ ), 



#200 = 0.156JV^ft _1 Mpc. 



(4) 



This relation was measured from the MaxBCG cluster 
sample, which ranged in redshift from 0.1 to 0.3 and had 
a median mass of about 1 x 10 14 h~ x M . The clusters 
we present in this work are mostly above this redshift 
range, and expected to have higher median mass (see 



V10, £4.2). Nonetheless, we adopt this relation to esti- 



mate cluster radii, and leave verification of the relation 
on an SZ selected sample such as this to future work. 

The entire richness procedure is repeated, now setting 
R = R200 instead of 1 h~ x Mpc, to arrive at an estimate 
of A^oo- A/200 is then used to estimate cluster mass using 
previously established empirical relations, which we now 
outline. 

3.3.1. Mass-Richness Scaling and Scatter 

By comparing to weak gravitati onal lensing masses , 
A/200 richness has been shown by |Reyes et al.| ([2008 ) 
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to scale with M2007 the cluster mass contained within 
a sphere that has an average mass density of 200 times 
the universal average, as 



M' 



200 



(1.42 ±0.08) 



200 



20 



1.16±0.09 



(5) 



Earlier, [Johnston et al. (2007) presented a similar, inde- 
pendent weliFTelisTn^stu^yT^nd their relation (normal- 
ization 0.88 ± 0.12, power law 1.28 ± 0.04) differs from 
Equation (pH) by about 30% in mass at 5 x 10 14 hr 1 M . 
We adopt Equation ([5| as our mass-richness scaling re- 
lation with overall uncertainty of 30%. 

Scatter in M200 at fixed A^oo = 40 as determined from 
weak lensing and X-ray cluster masses is <Jm\n — 45% ± 
20% (95% CL, [Ro"zo et al.(2QQ9| >. 

If richness quantities scale as 



200 



N, 



l/a 



200 



gal 



1/0.56 



gal 



R 



200 •> 



(6) 



where a is the N200-M200 power law and j3 is the AT ga i- 
A^oo power law, then a = 0.86 and f3 = 0.65, using re- 
sults of Equations @ and Scatter of 45% (60%) 
in mass translates to scatter of 39% (52%) in A/200 •> 
25% (34%) in JV ga i, and 15% (20%) in R 2m . 

As with the Agai-^200 relation (Equation Q), these 
7V~2oo~^200 mass relations were determined from weak 
lensing measurements of MaxBCG clusters in SDSS, a 
sample that was deemed to be complete at redshifts 
0.1 < z < 0.3, with median mass of approximately 
1 x 1O 14 /i _1 M . Our SZ selected cluster sample has 
a significantly broader redshift distribution, as well as a 
median mass about five times larger. Red galaxies are 
known to be biased tracers of dark matter, and the bias 
is a function of mass, redshift, and radius. Given the 
small size of our sample and the large intrinsic uncer- 
tainties inherent to richness techniques, we assume our 
data are insensitive to deviations from these fiducial re- 
lations, and we leave the study of the evolution of red 
galaxy populations with larger SZ cluster samples to fu- 
ture work. 

3.3.2. Statistical Richness Uncertainties 

We estimated sta tistical unce rtainties on N ga \ and A^oo 
by bootstrapping (Efron 1979) the entire richness pro- 
cedure thousands ol times! m particular, we assumed 
that the statistically limiting step in our procedure was 
the luminosity function fitting to the binned magnitude 
data. Therefore, the start of the bootstrap process was 
a random resampling of the magnitude bins themselves. 
For each realization, we re-fitted the luminosity function, 
integrated to obtain AT gal , and estimated i?200- We inde- 
pendently bootstrapped the LF fitting for A200 estima- 
tion as well. Parameters were assigned as the biweight 
means 
sample 1 



( |Beers et al.|19 90) of bootstrap distributions, and 
3 uncertainties as the biweight standard deviations. 



3.4. Spectroscopic Redshifts 

For each galaxy spectrum, the redshift was found by 
cross-correlating with the FA BTEMP97 template, using 
the RVSAO package in IRAF ( |Kurtz fc Mink|1998[ ). The 
validity of the cross-correlation redshift was checked by 
visual inspection and judged by the presence of visible 
absorption (and in a few cases, emission) lines. Redshift 




0.4 0.6 0.8 1.0 
Spectroscopic Redshift 

Fig. 5. — Spectroscopic versus red-sequence redshifts. We have 
applied an empirical linear correction to the red-sequence model 
colors using this sample, and this plot shows the result of red- 
sequence redshift measurements after the model correction. The 
best fit line to the uncorrected redshifts is the dashed line, shown 
for comparison. Typical RMS redshift scatter is about 2% in 
a z /(l + z). Redshift estimates for the entire sample are presented 
in Table [l| 



uncertainties were estimated as two times those given by 
RVSAO fo r the BCG redshift s, or the biweight interval 
estimator ( |Beers et al.||1990| ) for the other cases. We 
discard non-galactic spectra as well as redshifts in strong 
disagreement with the ensemble average or our prior pho- 
tometric redshift estimate. The redshift adopted for each 
cluster is the median redshift of the galaxies passing all 
cuts. 

4. RESULTS 

4.1. Redshifts 

We list redshift results in Table [I] Figure [5] is a plot of 
red sequence redshifts versus spectroscopic redshifts for 
clusters on which we have both measurements. The line 
we fitted to uncorrected redshifts, which were used t o em - 
pirically calibrate the red-sequence model colors (j |3.1| ), 
is the dashed line in the figure. After the correction, 
we verify that the red sequence redshifts are unbiased 
to within the uncertainties. Root-mean-square scatter in 
redshift per cluster is <j z =0.02(1 + 2), with maximum 
absolute deviation of 0.05(1 + z). We therefore generi- 
cally assign random uncertainties of 2% to all redshifts in 
the range we have directly tested, 0.15 < z < 1. For red- 
shifts estimated to be > 1, we assign 5% uncertainties in 
order to effectively rule out redshifts less than 1. Finally, 
for clusters at z < 1 where our photometric uncertain- 
ties are large, we also generically increase the uncertainty 
estimates. 

4.2. Richness 

In Figure [6] we compare our richness-derived masses 
to those presented in V10, which were estimated from 
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TABLE 1 



Cluster Redshift Data 


Cluster Name 


^rs 


^spec 




Imaging Coverage 








BCS? 


Magellan? 


SPT-CL 


J0509-5342 


U.4 1 {o) 


n a a \ 
U.4ozo^4 ) 


D 


Y 


Y 


SPT-CL 


J0511-5154 


U. i 4(^o ) 






N 


Y 


SPT-CL 


J0516-5430 


n o^t'Vi 

U.ZO^O J 


u.zyoz 


o 
o 


Y 


Y 


SPT-CL 


J0521-5104 


U. ( Z{6) 






Y 


N 


SPT-CL 


J0528-5259 


n 7K(a\ 


U. ( D4o^O ) 


o 
Z 


Y 


Y 


SPT-CL 


J0533-5005 


n &i.{a\ 

U.O<J^4 j 


n ssi n?Q^ 
u.ooiu^y ) 


/I 


N 


Y 


SPT-CL 


J0539-5744 


0.77(4) 






N 


Y 


SPT-CL 


J0546-5345 


1 1 C?1 1 \ 
1. lO^llJ 






Y 


N 


SPT-CL 


J0551-5709 


n /ii 
U.41(^o J 


U.4zoU(^lU J 


cr 
O 


N 


Y 


SPT-CL 


J0559-5249 


U.DO^O J 


U.Dllz^o J 


r 



N 


Y 


SPT-CL 


J2259-5617 


0.16(2) 


0.1528 


1 


N 


Y 


SPT-CL 


J2300-5331 


0.29(3) 






N 


Y 


SPT-CL 


J2301-5546 


0.78(9) 






N 


Y 


SPT-CL 


J2331-5051 


0.55(3) 


0.5707(5) 


8 


N 


Y 


SPT-CL 


J2332-5358 


0.32(3) 






Y 


Y 


SPT-CL 


J2337-5942 


0.77(4) 


0.7814(5) 


2 


N 


Y 


SPT-CL 


J2341-5119 


1.03(4) 


0.9983(5) 


1 


N 


Y 


SPT-CL 


J2342-5411 


1.08(10) 






Y 


N 


SPT-CL 


J2332-5521 






Y 


Y 


SPT-CL 


J2355-5056 


0.35(3) 






N 


Y 


SPT-CL 


J2359-5009 


0.76(4) 






N 


Y 


SPT-CL 


J0000-5748 


0.74(9) 






N 


Y 



Note. — See |4.3| for notes on individual clusters, including those which have 
been identified in other works. 



SPT millimeter-wavelength data. We label V10 masses 
as M(£), because they are calculated from the SZ signal- 
to-noise ratio £. The purpose of the comparison is to 
assess the level of correlation and explore whether the 
power law and normalization of the scaling laws are con- 
sistent with previous work. We additionally estimate the 
normalization of iV" ga i, which we treat as an empirical 
mass observable in its own right rather than only as a 
measurement intermediate to A^oo- 

Richness and M(£) are plotted against one another in 
Figure [6j We have used Equation (|5| to render the far 
right-hand mass axis, with the 30%overall uncertainty 
in this relation denoted with the heavy error bar. We 
do not display a mass axis in the N g3 \ panel of the figure 
because there have been no previous measurements of an 
7V ga i-M2oo scaling relation. 

We perform the analysis only considering clusters 
whose radius (R200 m the case of A200? lft -1 Mpc in the 
case of iVg a i) falls fully within the observed field of view. 
For Magellan IMACS imaging this limited radii to < 6', 
because we placed the SPT center in the middle of one 
of the eight chips, at a distance of about six arcminutes 
from the field edge. Spatial incompleteness of this kind 
does not affect the BCS data. 



4.2.1. Agai Scaling 

i?2oo f° r massive clusters is approximately 1.5 ft -1 Mpc, 
which subtends > 6' at redshift s z < 0.3. We find it 
useful to count red galaxy overdensities in a smaller re- 
gion, corresponding to 1 h~ x Mpc in co- moving coordi- 
nates, which subtends > 6' at redshifts z < 0.17. This 
is a better match to the observed typical angular size 
of the SZ signal. In addition, the surface density of red 
galaxies in the smaller aperture is greater than for i?200, 
whose larger area effectively dilutes the signal. For these 
reasons we explor e A^ ga i as an empirical mass proxy. 

As argued in § 3.3.1 , Af ga i should scale with mass as 



Mijoo 6 - We fix this slope and measure 



A^ ga i = (52 ± 3 ± 9) 



200 



5 x 10 14 HM 



0.56 



(7) 



If zero intrinsic scatter is assumed, the reduced chi- 
square of the fit is xt ~ 3, whereas 20% in N ga \ (cor- 
responding to ~ 35% in mass) produces xt = 1- The 
first error given in Equation (I7|) is random-only, which 
includes the intrinsic mass-richness scatter, and the sec- 
ond error is the overall systematic uncertainty of 30%. 
This fit is shown as the solid black line in the left-hand 
panel of Figure [6j 

Letting the slope be free, we measure a normalization 
consistent with Equation (ITl and slope 0.42 ± 0.50. This 
is the dashed line in the lert-hand panel of Figure [6] In 
the figure, the inner, dark shaded area denotes the 68% 
confidence region for the best two-parameter fit assum- 
ing zero intrinsic scatter, while the outer, light shaded 
area is the 68% confidence region for best two-parameter 
fit assuming 35% intrinsic mass scatter, which is a sig- 
nificantly better match to the data. We estimate the 
uncertainty on intrinsic scatter by varying it until the 
chi-square doubles, resulting in relative la uncertainty 
at the - 50% level. 



4.2.2. N 2 oo Scaling 

We perform the same analysis on A' 2 oo- 
mass-mass slope to unity, we measure 



Fixing the 



^200 = (57 ± 4 ± 15) 



AL 



200 



5 x K^ft-iM© 



0.86 



(8) 



Again, 35% intrinsic mass scatter, corresponding to 30% 
scatter in A/200 > gives xt = 1- This fit is shown as the 
black line in right-hand panel of Figure [6| and can be 
compared to the one-to-one mass line, in red. 

Letting the slope be free, the measured normalization 
consistent with Equation ([8| and slope 0.62 ±0.73. This 
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Fig. 6. — Optical richness versus total cluster mass estimated from the SPT dat a. This plot shows that the richness correlates highly with 
this millimeter- wavelength mass-observable (taken from Vander linde et al.|2010j ), and together the data agree with previously established 
scaling relations to within the uncertainties. The solid red line m the right-hand panel is the one-to-one mass relation, and the solid black 
lines in both panels are best-fit relations when fixing the slope to previously measured values. The dashed lines are best-fit relations leaving 
both slope and intercept free. The dark, inner shaded areas denote the 68% confidence regions assuming zero intrinsic richness-mass scatter 
for the two parameter regression. The light, outer shaded areas denote 68% confidence regions assuming 35% scatter in mass, which is a 
better match to the data. The heavy error bar on the far right-hand side indicates the nominal overall uncertainty of 30% in the A^200~^200 
relation. 



is shown with the dashed black line in Figure [6j As be- 
fore, the dark, inner shaded region is the 68% confidence 
region for the best two-parameter fit assuming zero in- 
trinsic scatter, while the outer, light shaded area is the 
68% confidence region for best two-parameter fit assum- 
ing 35% intrinsic mass scatter, which is a significantly 
better match to the data. Relative uncertainty on the 
intrinsic scatter is comparable to that quoted above for 

All richness results aregiven in Table [2] In the table we 
have adopted Equation (pi to estimate masses from A/200 ? 
and for N g3 \ masses we nave used Equation ([7]). System- 
atic uncertainties are taken to be the quadrat ure su m of 
the nominal 45% intrinsic scatter in mass (j |3.3.1 ) and 
the 30% overall uncertainty in the richness-mass scaling 
relation. Further work on a larger sample of clusters se- 
lected with similar criteria as these is needed to reduce 
statistical uncertainties and measure the scatter directly. 
We discuss the implications of our richness measurements 
in© 

4.3. Notable Clusters 

In this section we describe notable information, if any, 
about the clusters. As we will point out, a subset of 
clusters also appear in the catalogs of |A89| [B ohringer 



overlapping clusters is too small to draw useful conclu- 
sions with high statistical confidence, we leave a fo rmal 
inte r-com parison of redshift and richness results of M09 
and IMlQl to future studies. 



SPT-CL J0509-5342- 
tified by[SQ9l 



SPT-CL JO 5 1 1-5154 
identified by |M10 



J051145-515430. 



This cluster was previously iden- 



— This cluster has recently been 
who assigned the name SCSO 



SPT-CL J0516-5430— This cluster was previously iden- 
tified by S09, where it was called by a different name, 
SPT-CL J0517-5430. The SPT name ascribed to this 
object in this work and in VI follow the recommenda- 
tions of the International Astronomical Union (IAU) , and 
should be adopted pe rmanently. This c luster also iden- 
tified as Abell S0520 dAbell et al.1|1989[ hereafter A89), 



and RXCJ0516.6-5430 ( IBohringer et al.|2004| ), the latter 
of which is the source 01 the spectroscopic redshift. [M1Q| 
detected this object and call it SCSO J051637-543001. 

S PT-CL J0521-5104— This cluster has been identified 
by |MlQl as SCSO J052113-510418. 



et al.l (|2QQ4l |SQ9l [M enantea u fc Hughes | ( |2QQ9["h"ere - 
after MH09), Menanteau et al. (2009, hereafter M09), 
and |Menanteau et al.| ( |2010| hereaFter |M10| ). OurTed- 
shifts agree to within <t7" ^~3% (1 + z) ol the photomet- 
ric redshifts presented in M09 ] and |M10[ except for the 
highest redshift cluster from |S09[ SPT-CL J0546-5345 
(see below). Because we have not presented exhaustive 
optical cluster-finding in the entire BCS survey in this 
work, instead having concentrated on fields in the direc- 
tion of the SZ detections, and because the number of 



SPT-CL J0528-5300— This cluster was previously iden- 
tified by S09, where it was c alled by the same name, and 
was also identified by [M10l as SCSO J052803-525945. 

SPT-CL J0539-5744— This cluster displays a possible 
strong gravitational lens arc. 

SPT-CL J 0546-5345— This cluster was previously iden- 
tified by S09, where it was called by a different name, 
SPT-CL J0547-5345. The SPT name ascribed to this 
object in this work and in VI follow the recommenda- 
tions of the IAU, and should be adopted permanently. 
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TABLE 2 

Cluster Richness Data 



Object Name 


^gai a 


M 200 (^gai) b 


#200 a 


A^200 a 


M 20 o(A^20o) b 






(lO^/i^M©) 


(/i _1 Mpc) 




(lO 14 ^- 1 M©) 


SP1-CL J0509-5342 


41(8) 


3.3 ± 2.0 ± 1.0 


1.46(17) 


51(13) 


4.3 ± 2.5 ± 1.3 


bFl-CL J0511-5154 


77(12) 


10.3 ± 5.9 ±3.1 


2.13(20) 




SP1-CL J0516-5430 








SP1-CL J0521-5104 


44(7) 


3.7 ± 2.2 ± 1.1 


1.52(15) 


57(13) 


4.8 ± 2.7 ± 1.4 


nnrri pj Tnroo conn 

br 1-UL J0528-5300 


44(9) 


3.0 ± 2.3 ±1.1 


1.52(lo) 


60(13) 


EC 1 _L O O _l_ 1 EC 

5.1 ± 2.8 ± 1.5 


SP1-CL J0533-5005 


28(10) 


1.7 ± 1.4 ± 0.5 


1.17(25) 


28(6) 


z.z ± 1.2 ± 0.6 


bPl-CL J0539-5744 


63(9) 


7.1 ± 4.0 ± 2.1 


1.88(16) 


69(14) 


1 100110 
o.l ± 3.3 ± 1.8 


br 1-UL J054o-5345 


66(7) 


7.7 ± 4.1 ± 2.3 


1 no/1 o\ 
1.93(12) 


80(31) 


7.1 ± 4.8 ± 2.1 


bFl-OL J0551-5709 


54(15) 


5.4 ± 3.8 ± l.b 


1 71 /on\ 

1.71(29) 






SPT-CL J0559-5249 


59(6) 


6.5 ± 3.4 ± 1.9 


1.82(11) 


70(8) 


6.2 ± 3.2 ± 1.8 


SPT-CL J2259-5617 












SPT-CL J2300-5331 


35(10) 


2.5 ±1.8 ±0.8 


1.33(22) 






SPT-CL J2301-5546 


35(8) 


2.5 ±1.6 ±0.8 


1.33(17) 


42(6) 


3.4 ±1.8 ±1.0 


SPT-CL J2331-5051 


73(5) 


9.4 ±4.8 ±2.8 


2.06(8) 






SPT-CL J2332-5358 


42(8) 


3.5 ±2.1 ±1.1 


1.49(16) 






SPT-CL J2337-5942 


53(8) 


5.2 ± 2.9 ± 1.6 


1.69(15) 


55(11) 


4.6 ±2.5 ±1.4 


SPT-CL J2341-5119 


39(7) 


3.1 ±1.9 ±0.9 


1.42(16) 


62(7) 


5.4 ±2.8 ±1.6 


SPT-CL J2342-5411 


43(19) 


3.7 ± 3.4 ± 1.1 


1.51(38) 


37(12) 


2.9 ± 1.8 ±0.9 


SPT-CL J2355-5056 


55(5) 


5.6 ± 3.0 ± 1.7 


1.73(10) 






SPT-CL J2359-5009 


47(16) 


4.3 ± 3.4 ± 1.3 


1.59(32) 


53(8) 


4.4 ±2.3 ±1.3 


SPT-CL J0000-5748 


45(6) 


4.0 ± 2.2 ± 1.2 


1.55(12) 


54(11) 


4.6 ±2.5 ±1.4 


Note. — 



Uncertainties given are statistical only. 

Uncertainties given are statistical and systematic, respectively. 



[SQ9] reported a photometric redshift of - 0.9 and M09 

' ' _i_n n& ' ' 

independently reported z p hoto = 0.88_ ; 04 . We do not 
detect a red-sequence overdensity near redshift 0.9, but 
a faint red-sequence peak is evident at z rs ~ 1.15. 

SPT-CL J0551-5709— Abell S0552 flA89| ) is in the fore- 
ground of this cluster. No redshift estimate exists for 
Abell S0552, but a strong red sequence at z rs = 0.09 is 
clearly visible in color- magnitude diagrams. The SPT 
cluster identified here, however, is measured at z = 0.42. 

SPT-CL J0559-5249 — Our red-sequence studies reveal 
two significant red galaxy over densities, one at the red- 
shift given in Table [T] and another at z ~ 0.4. On inspec- 
tion of the spatial distribution of galaxies, we attribute 
the SZ signal to the higher redshift system. 

SPT-CL J 2259-5611— We identify this cluster with 
Abell 3950, and recover a spectroscopic redshift from 
archival data on what we identify as the BCG. The SPT 
SZ detection coordinate lies nearly exactly on the line 
joining the two Abell 3950 coordinates given in the lit- 
er ature, at a projected d istance of 71" from that quoted 
by |Arp fc Madore| ( |1996T >, and 208" from that quoted by 
A89. The redshilt ol Abell 3950 has not been previously 
measured, but we have identified the BCG in 2MASS 
as 2MASX J23000108-5617061, which the 6dF Galaxy 
Surve y measured to be at z spec = 0.152787 (Jones et al. 
|2005[ ). This galaxy lies 18" from the SPT coordinate. 

SPT-CL J2300-5331— We identify this with Abell S1079 
(A89). No previous redshift estimates exist for this clus- 
ter. 

SPT-CL J2331-5051— This cluster exhibits a giant grav- 
itational lens arc and a well separated secondary cluster 
structure in both the optical and SZ data. This is among 
the most interesting of the clusters presented here, and 



is the subject of a dedicated study (High et al., in prepa- 
ration) . 



SPT-CL 



'-5358 — This cluster was recently identi- 



fied by |MlQ| as SCSO J233227-535827. 



SPT-CL J2343-5521 — No red-sequence cluster appears 
in BCS imaging, whose 50% completeness depth in the 
i band is 23.5 mag, corresponding to m* at z ~ 1.2. 
Either there exists a cluster at higher redshift, or this is 
a false SZ detection. V10 show that the false detection 
rate is only 7% for clusters at equivalent SZ significance. 
However, the SZ profile radius is significantly larger than 
any other cluster, consistent with a cosmic microwave 
background fluctuation. Preliminary follow-up of this 
cluster at other wavelengths suggest this is indeed a false 
detection. 

5. DISCUSSION OF SYSTEMATIC EFFECTS 

The greatly different criteria with which our SZ clus- 
ters were selected as compared to the MaxBCG sample 
could give rise to differences in measured properties. One 
important effect is the evolution of the mean and scat- 
ter of red cluster-galaxy colors with redshift and mass. 
While elliptical (E) galaxies are highly homogeneous in 
the range of redshifts our sample represents, < z < 1. 2 
flMenci et aL|2008| |Lidman et aT|2008} |Mei et al.|[2 009 ) , 
the mean and scatter ol SO galaxy colors are known to 
evolve with redshift and d ensity of the environment (e.g., 

The environment also causes 
from 



van Dokkum et al 
color evolution 



1998) 
witrT 



distance from the cluster center. 
Our method, and indeed that of MaxBCG, does not ex- 
plicitly use morphological selection criteria, so we must 
assume E and SO populations contribute to our richness 
estimates. Understanding color-selection effects at all 
redshifts for a large sample of clusters would require so- 
phisticated simulations or full photometric redshift mea- 
surements that are beyond the scope of this work. 
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Another important effect is the evolution of the abun- 
dance of early-type cluster galaxies as a function of lumi- 
nosity. Clusters have been observed to accumulate faint 
galaxies at a greater rate than bright ones over time, 
manifesting as evolution in the slope of the luminosity 
function's faint end, a, as a function of redshift ( |Rud-| 
nick et al.||2009|). Despite our long redshift baseline, we 



expect our measurements here to be largely unaffected by 
behavior at the faint end, as we integrate the LF down 
to a relatively shallow magnitude of ra* + 1. Indeed, 
tests where we vary a have not significantly affected our 
results. 

In addition, because our richness measurement is not 
exactly that used with maxBCG (from which the i^200~ 
TVgai relationship that we use here is taken) our i? 2 oo 
estimates are likely to be somew hat overestimated (cf. 



discussion in Hansen et al. 



2009). Using a larger-than- 



ideal aperture lor counting red galaxies may contribute 
to some of the scatter in the richness-mass correlation 
that we observe here. 

Red galaxy counting has nonetheless proven to be a 
simple and accessible way to estimate the total mass in 
clusters and groups. If this technique can be accurately 
extended to the very wide range of redshifts that SPT 
SZ-selected clusters span, then modeling and measuring 
evolutionary effects in SZ clusters will be useful to ob- 
taining constraints on masses of very large cluster sam- 
ples. 

6. CONCLUSIONS 

We have observed clusters from the 2008 South Pole 
Telescope SZ survey at optical and near-infrared wave- 
lengths. We estimate redshifts and richness with red- 
sequence techniques, and we obtain spectroscopic red- 
shifts for a subsample of the clusters. 

Our red-sequence-derived redshifts exhibit 2% RMS 
scatter in cr z /(l + z) in the subsample with spectroscopic 
overlap, over the redshift range 0.15 < z < 1.0. Our anal- 
ysis provides no evidence that the SZ selected sample 
from SPT follows different scaling relations than those 
followed by SDSS optically selected clusters. 

The clusters presented in this paper comprise the 
largest sample of galaxy clusters discovered with the SZ 
and demonstrates that current SZ surveys can detect 
many high-mass galaxy clusters across a wide range of 
redshifts. Precise dark energy constraints from these sur- 
veys require the cluster redshifts, masses, and selection 
function to be known. The SZ effect contains no redshift 



information, and coordinated observations at optical and 
infrared wavelengths are an efficient means of providing 
this, especially for large cluster samples. Optical cluster 
identification is also potentially useful for understanding 
SZ selection functions. In the future, there will proba- 
bly be greater coverage of SZ clusters by OIR imaging 
than any other wavelength or observing method other 
than the millimeter-wavelength itself. Optical redshifts 
are essential for constraining cosmology with SZ surveys, 
and the same data may also be brought to bear on the 
mass and cluster selection problems. 
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APPENDIX 
GALLERY OF CLUSTERS 



Figures [7[j28| are false-color optical images of the clusters, together with the SZ detection significance maps. In all 
images, northis up, east is left. 
The SZ-only insets subtend 8 arcminutes on a side. The mapping between color and SZ significance £ is the same in 



all SZ thumbnails. The highest significance cluster, Figure 22, spans the largest color range , which is £ G [—5.57, 14.94]. 
The peak value in each thumbnail is equal to the quoted SZ detection significance in |V10| Contours are given at steps 
of A(£) = 2, and the corresponding values are labeled. Contours are dashed where £ is negative, and solid where £ is 
nonnegative. 

The optical images have the same contours as their corresponding SZ thumbnail overlain, and subtend 1.5 Mpc in 
the cluster rest frame. We have mapped the zrg passbands to the RGB channels, respectively, so that the apparent 
redness of red-sequence cluster galaxies approximately increases with redshift. 
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Fig. 8. SPT-CL J0511-5154 




Fig. 10. 



— SPT-CL J0521-5104 
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Fig. 12. SPT-CL J0533-5005 
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Fig. 13. SPT-CL J0539-5744 



1 arcmin 




Fig. 14. 



SPT-CL J0546-5345 
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Fig. 15. SPT-CL J0551-5709 




Fig. 16. SPT-CL J0559-5249 




Fig. 18. SPT-CL J2300-5331 
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Fig. 19. SPT-CL J2301-5546 




Fig. 20. SPT-CL J2331-5051 




Fig. 22. SPT-CL J2337-5942 
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Fig. 23. SPT-CL J2341-5119 




Fig. 24. SPT-CL J2342-5411 
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Fig. 28. SPT-CL J0000-5748 



